Exact Kohn-Sham exchange kernel for insulators and its long-wavelength behavior 
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We present an exact expression for the frequency-dependent Kohn-Sham exact-exchange (EXX) 
kernel for periodic insulators, which can be employed for the calculation of electronic response prop- 
erties within time-dependent (TD) density-functional theory. It is shown that the EXX kernel has 
a long-wavelength divergence behavior of the exact full exchange-correlation kernel and thus recti- 
fies one serious shortcoming of the adiabatic local-density approximation and generalized-gradient 
approximations kernels. A comparison between the TDEXX and the GW-approximation-Bethe- 
Salpeter-equation approach is also made. 
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I. INTRODUCTION 

Timc-dcncndcnt density-functional theory 

(TDDFT)Era is an attractive first-principles formal- 
ism for the calculation of electronic response propertie; 
Thanks torits simplicity, applications to quantum wellsl 
and atomsu have appeared long before its formal justi- 
ficationja and currently its use is ubiquitous in a wide 
range of ah initio computations. A crucial ingredient of 
successful TDDFT applications is the approximation to 
the dynamic exchange-correlation kernel, 



F xc (r,r';t-t') 



faxc(r;Q 
5n(r';t') ' 



(1) 



which together with the Hartree kernelH i%(r, r') = 
Wc(r, r') = 1/ 1 r — r'| completely determines the two- 
particle interaction effects. For F xc , the adiabatic local- 
density approximation (LDA) kernel, 



F 



LDA 



(r,r';t 



dn(r) 



(2) 



has been almost exclusively adopted in practical calcu- 
lations. However, the scope of the LDA kernel has been 
rather limited for infinite periodic solids due to its de- 
ficiencies, and in particular its incorrect non-divergent 
long-wavelength behavior for insulators has-been empha- 
sized as a primary defect in recent years. Qtm This short- 
coming of the LDA kernel shows up, e.g., in its incapa- 
bility of describing excitonic effects in absorption spectra 
of solids. The semilocal generalized gradient approxima- 
tions (GGA) kernel does not improve over the LDA one 
in this case, and the task of developing a more accu- 
rate approximate exchange-correlation kernel remains as 
a challenging task for the TDDFT study of solids. 

Indeed deficiencies of the LDA and GGA appear al- 
ready at the level of the static exchange-correlation en- 
ergy functional E xc [n] and the exchange-correlation po- 
tential v xc (r) = 8E xc [n)/5n(r). For instance, the LDA 
and GGA E xc [n] inherently fails to describe the quasi- 
two-dimensional electron gas due to their (semi)local na- 
ture.Ej For v xc and the corresponding Kohn-Sham (KS) 
eigenvalues, the LDA and GGA v xc incorrectly decays 



exponentially rather than as —1/r for localized systems, 
and consequently their highest occupied orbital energies 
are too high and unoccupied orbital energies do not ex- 
hibit Rydberg series. For solids, the LDA band gaps are 
too small and this behavior is again not corrected by 
the GGA. In fact, even the exact KS gap does not equal 
the experimental band gap but differs by the disconti- 
nuity of the exchange-correlation potential.til However, 
recent theoretical and numerical studies suggested that 
the KS equation for N electrons corresponds to the Dyson 
equation where tiie reference ground state is chosen with 
N — 1 electrons,E3 and accordingly unoccupied orbitals 
in the KS calculations should give a good description of 
excitations of the ^-electron system with fixed particle 
number. This is in accordance with—the perturbation 
theory along the adiabatic connection^ which finds that 
differences of KS eigenvalues represent the leading term 
in the expansion of excitation energies. 

In this regard, recent development of the KS exact- 
exchange (EXX) method, which treats the exchange- 
correlation energy functional exactly in leading order 
in the electron-electron interaction, provides an inter- 
esting opportunity. Self-interaction-free, nonlocal EXX 
schemes give not only realistic exchange potentials and 
KS eigenvalue spectra for moleculesliJ'EjIlj'EJEa but also 
band structures of ..semiconductors in good agreement 
with experiments. EjE3 We have recently shown that the 
EXX orbitals and eigenvalues at the one-particle level 
without any previously applied post-DFT modification 
such as the quasiparticle shiftEil indeed give a very good 
description of the absorption spectrum of semiconductors 
with the exception of excitonic features resulting from 
two- particle interactions, and argued that it is another 
evidence, of the above-described picture of "KS quasipar- 
ticles" £3 

In view of the encouraging performance of the EXX 
method, we present in this work an exact expression of 
the EXX kernel F® for periodic insulators which can 
be employed for calculations of electronic linear response 
properties within TDDFT. It will be shown that the EXX 
kernel, unlike the LDA and GGA kernels, exhibits a long- 
wavelength behavior of the exact F xc which is particu- 
larly important for the study of electronic excitations in 
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infinite solids. This behavior of the EXX kernel has been 
previously claimed by Goshez et al. based on a plausi- 
bility argument ,@ and here we explicitly prove this using 
our exact formula. 



II. LONG- WAVELENGTH BEHAVIOR OF THE 
EXCHANGE-CORRELATION KERNEL FOR 
INSULATORS 

We first establish the definitions and notations of vari- 
ous quantities of interest and derive the long-wavelength 
behavior of the exact F xc .eI The full linear density re- 
sponse matrix x describes the response of the first-order 
(number) density change Sn for the given bare dynamic 
perturbation cH> ox t, 



Sn(G, q; u) = £ X (G, G', q, u) 6v ext (G', q; u). (3) 

G 



We chose to work in the reciprocal space and the fre- 
quency domain, e.g., x(G,G',q;w) is a matrix in the 
reciprocal-space lattice vectors G and G' for the given 
wave- vector q and frequency to. Analysis of the q — > be- 
havior of the "head" (G = G' = 0), "wing" (G = and 
G' ^ or vice versa), and "body" (G ^ and G' ^ 0) 
elements of x and other related matrices appearing below 
is an important discussion point throughout the paper. 
From now on, we will adopt the matrix notation and G 
dependence will be assumed unless explicitly stated oth- 
erwise. Within TDDFT, 8n is expressed in terms of the 
dynamic linear response matrix xo and the first-order 
change of the effective KS potential Svks, 



5n(q;uj) = xo(q; w) <5v K s(q; w ), 



(4) 



where Svks is composed of the external perturbation 
<5i> ex t and the resulting change in the Hartree potential 
8vn and the exchange-correlation potential 5i> xc , 

fe K s(q; w) = <5v ox t(q; w) + <5v H (q; w) + <5w xc (q; ui) 

= <5w cxt (q; uS) + [F H (q) + F xc (q;w)] <5n(q;w), 

(5) 

with F H (G,G' q) = (5 G , G '47r/|q + G| 2 . Then, from Eqs. 
(|), (§), and (g), one obtains 

Xo x (q; w) = x _1 (q; w) + F H (q) + Fcc(q; w), (6) 

which shows that x is completely determined once xo and 
F xc are given. 

For further consideration of the q — > behavior of Fxc, 
it is convenient |-to-,introduce the "proper" part of x, Xj 
defined throughBE3 

5n(q\u) =x(q;w)^ TC (q;w), (7) 

where <5i>tc is the change of the test-charge potential, 

<5v T c (q; w) = <5w ext (q; w) + F H (q) 5n(q; w), (8) 



and thus relates with the full response matrix x as 

X -1 (qjw) =X _1 (q;^)+^H(q), (9) 
or with the KS response matrix xo as 

Xo 1 (q;w) = X" 1 (q;w)+F xc (q;u;). (10) 

The linear response matrix of the non-interacting KS sys- 
tem xo and the proper part of that of the real interacting 
system % are , known to have the following similar q — > 
behaviorEjE.il (assuming that both the KS system and the 
real system are insulatingu) : 



Xo 



q 2 X°o° QXa 1 

rf xh 1 



X 



q 2~00 q ~01 

qX W X 11 



(11) 



where we used the notation that £ 00 , £ 01 / 10 , and £ n de- 
note the head of a matrix £ divided by q 2 , its wings di- 
vided by q, and its body, respectively. The quantities 
Xo°, xl°, Xo 1 , and x" as well as x 00 , X 10 , X°\ and X U 
all consist of a leading (^-independent term and contribu- 
tions of higher order in q which vanish in the limit q — » 0. 
Then, from Eqs. (10) and (|TT|) , one can deduce that the 
head ancLwings of F xc have the following divergent q — > 
behaviors (assuming that there exists no fortuitous can- 
cellation between Xq 1 and xT 1 ), 



F™/q 2 F£/q 

F w /q F 11 
xc / y xc 



(12) 



In Eq.@, F° c °, F™, F°£, and F 11 again contain a lead- 
ing ^-independent term and contributions of higher order 
in q which vanish for q — > 0. The head and wings of adi- 
abatic LDA and GGA kernels, on the other hand, are 
independent of q and thus are incorrectly non-divergent 
for q — > 0: 



pLDA/GGA 



p LDA/GGA,00 j-,LDA/GGA,01 

xc xc 

pLDA/GGA, 10 pLDA/GGA, 11 

xc ^xc 



(13) 



This defect is a serious problem not only from a theoreti- 
cal viewpoint but also for practical purposes because the 
head and wings of the exchange-correlation kernel can 
affect Jjhe the macroscopic dielectric function in leading 
order .0 For example, it has been recently shown that they 
play a crucial role for the proper treatment, of excitonic 
effects in the calculation of optical spectra.S 



III. EXACT-EXCHANGE KERNEL AND ITS 
LONG- WAVELENGTH BEHAVIOR FOR 
INSULATORS 

Deficiencies of the LDA and GGA kernels discussed 
above represent a major problem from the theoretical and 
calculational point of view which could not be overcome 
so far. To ameliorate the situation we propose to adopt 
the EXX kernel. An exact expression of the EXX kernel 
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has been previously derived by one of us-.for localized 
systems for the case of real-valued orbitals.Ej For periodic 
solids, we need to generalize this expression to complex 
orbitals and have to consider the dependence on wave 
vectors k and q. This leads to 

F™(G,G,q;u;) = X^G.Gi.qjw) 

Gi,G 2 (14) 

x H X (G U G 2 , q; w) x " 1 (G 2 , G', q; «), 



where the EXX kernel "core" H K is composed of the 
following contributions (We assume that £u ext (q;w) 
and other quantities have the time-dependence e~ luJt e St , 
where 6 — > + is a convergence factor.): 



(ak|e- l ( c i+ G )- r |sk + q)(sk + q; bk'\w c \tk' + q; ak)(tk' + q|e*( q+G ')' r |&k') 
s k , k , - (e a k - esk+q + to + i£)(e&k' - e<k'+ q + oj + iS) 

(sk|e-*( q+G )- r |ak + q)(ak + q; tk'\w c \bk' + q; sk)(6k' + q|e l ( c ' +G ')- r |ik / ) 



(15) 



(e a k - e s k+ q - w - %S)(ebV - e*k'+q -u—iS) 



fl2(G,G' >q;W ) S -ix;E[ 



(ak|e-' i ( cl+G )- r |sk + q)(sk + q; ik'|w c |6k' + q; ak)(6k' + q|e^ q+G ')- r |ik') 



ask btk.' 



+ 



(e Q k - e s k+q + w + i^)(e&k' - e*k' 



a; — i(5) 



(sk|e-' i ( c '+ G > r |ak + q)(ak + q; &k'|uJ c |ik' + q; sk){tk' + q|e^ q + G ')- r |&k') 



(16) 



(e a k - e s k+ q - w — iS)(e bk > - e ik '+ q + w + iS) 



nun n> \- 2 \^ r (ak| e -^q+ G )- r lsk + qX&kj^ - v x \ak)(sk + q|e'(q+ G ')->k) 
S2 ^ L e„k - e.9k+n + w + ifl eik - e sk +n + w + io 



afrsk 



(e a k - e.sk+ q + w + «<5)(e&k - e s k+ q + lu + iS) 
(sk|e-' i ( c i+ G )- r |ak + q)(ak + q|S K - w^Mc + q)(6k + q|e 4 ( q+G ')' r |sk) - 



(e a k - Esk+q —oj — iS) (etk - e s k+q — u — iS) J 
2 ^ r (ak|e- l ( q+G )- r |sk + q)(sk + ql^ - v x \tk + q)(tk + qle^+^^ak) 
fi ^ L (e ak - e s k+ q + w + iS) (e ak - e tk+q + to + iS) 



(sk|e- i ( cl + G )- r |ak + q){tk\Y, x - v x \sk)(ak + q\e l ^ +G '> r \tk) 



(e a k - e s k+ q - us — iS) (e ak - e t k+ q — a — id) 



4 



and 



tf x 4 (G,G',q;c)^]r[ 



(bkle-^+^-'Isk + q)(sk + q\t x - v x \ak + q)(ak + q|e l ( q+G ')- r |6k) 
(e&k - e s k+q + lu + «<5)(e Q k - e s k+ q ) 

(ftkle-^^^-'lak + q)(ak + q|S x - d x \sk + q)(sk + q|e^ q+G > r |&k) 

(eak - £sk+ q )(efck - e s k+q + uj + i5) 
(sk|e- i ( q+G )- r |6k + q)(ak\t x - t> x |sk)(6k + q|e i ( q+G '> r |ok) 



in 



asik 



(Cbk - 


Csk+q — W 


- «<5)(£ak - £sk+q) 




(ak|e- 4 ( c '+ G )- r |6k- 


f q)(sk|S x 


-u E |ak)(6k + q|e 4 ( q+G '' 


• r |sk) 


(fiak - 


£ sk+q)(ebk 


- £sk+q —LJ—i5) 




(ak|e- 4 ( c '+ G )- r |ik- 


fq)(sk|E x 


-^|ak)(tk + q|e i ( q+G ') 




(e a k — 


£tk+q + W - 


f i£)(e a k - e s k+ q ) 




(sk|e- i ( q+G )- r |tk- 


hq)(ak|g a; 


-?5 x |sk}(tk + q|e i ( q + G ') 


r |«k) 



(e a k - e s k+q)(e a k - etk+q + lu + iS) 
(tkle-^i+^-'lak + q)(ak + q|E x - ^jsk + q)(sk + q|e l ( q+G ')- r |ik) 

(e a k — e tk+q ~ W — 4(5) (e a k — e sk+q) 

{ik|e-^ q+G >- r |sk + q)(sk + q|S x - S^ak + q)(ak + q|e 4 ( q+G ')- r |ik)- 



(18) 



In Eqs.@-(p 
(ik + q; jk'\w c 



( e ak ~ e sk+q)(eak — e tk+q ~ LO — iS) 

), 2 is the spin factor, is the crystal volume, {a, 6} are valence bands, {s,i} are conduction bands, 
Ik! + q; mk) are four-index integrals defined as 



{ik + q; jk'\w c \lk' + q; mk) = dr dr 



, ^k +q (r)^V(r')^k' +q (r)0 mk (r') 



4vr ^ _^k+_q|e^ G +^- k ')- r |Zk' + q)(jk'| e -' i ( G + k - k ' 



|r7ik) 



(19) 



|G + k-k' 



S K is a nonlocal orbital-dependent exchange operator of the form of the Hartree-Fock exchange operator but con- 
structed with the KS orbitals 6 a , 



(ik + q|S x |jk + q) = - JdrJ dr'# k+q (r) ^ 

nXe' 



</w(r)^ k ,(r') 



»jk+ql 



= -^^ k+q;ak 'i^ c, i ak ' ; - 7 ' k + q ) 



a k' 



4tt ^ (zk + q|e^ G+k ~ k,+q) - r |ak , )(ak , |e- i ( G+k - k, +q)- r '|jk + q) 



ak'G 



IG + k-k' 



r 



(20) 



and v x is generated by the local multiplicative EXX KS 
potential v x (r). 



Compared with the LDA (or GGA) kernel which 
is (semi)local in real space and frequency-independent 
[Eq. (0)], which results in a reciprocal- representation 



independent of q and u>, F, 

LDA/GGA, 



LDA/GGA 



(G,G',q;o;) 



F 5 £ 1 " l/l "" 1 (G-G , ) ) F® is fully nonlocal in real space 
and depends explicitly on the frequency. We now show 
that F^ xx has a q — > behavior as the exact F xc . By 
expanding orbitals ^k+q in terms of the orbitals em- 



ploying perturbation theory. 



Ojk+q 



>ik- 



q ■ (tkjpjjk) 

£ik — £jk 



(21) 



we express various matrix elements of Eqs.(|l5[) — (|Tq) in 
power series in q and keep only the leading non- vanishing 
terms. Then, one can first observe that (ik; jk' \lk' ; mk) 
and (ik\T, x — v x \jk) are the leading order terms in q of 
(ik + q; j'k'|ik' + q;mk) and (ik + c\\Y, x — v x \jk + q) and 
that consequently the q-dependence can be ignored for 
q — > in the inner matrix elements of Eqs. (|l5|) 
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One might notice that the inner matrix elements (iO + 
q; jO\iO + q; jO) in contain a singular contribution, 
the term with G = in Eq. (|l9]). However, the same 
singularities with the opposite sign arise in the matrix 
elements (ik + qlE^ljk + q) of the first two contributions 
of H x , the terms with G = and k' = k + q in Eq. 
©. So the Coulomb singularities in the inner matrix 
elements of H x and H x exactly cancel. 

Unlike in the case of inner matrix elements, for the 
other outer matrix elements (ik|e~ l< - q+G - l ' r | jk + q) or 
(ik + q|e~^ q+G • ) r |jk) with i ^ j, q-dependent contrihai- 
tions appear in leading order in q for G = or G' = OcJ, 
e.g., 



lim(ik|e~ 

<7—0 



iq-r 



lit- 



<ik|p|jk> 



(22) 



In H x , however, matrix elements (ik|e l ( q+G )' r |zk + q) 
and (ik + q|e _I< - q+G )' r |ik) are present, for which a lead- 
ing order term independent of q occurs for G = or 
G' = 0. However, contributions of such type in the first 
sum of are cancelled by corresponding contributions 
in the second sum. 

Due to the cancellations of singularities, H x itself is 
well-defined, and the q — > behavior of H x can be de- 
duced as 



H 



(23) 



with H x °, H*°, H® 1 , and H^ 1 containing a leading order 
term independent of q. Consequently, using Eqs. (11), 
(fj), and @, we conclude that F® xx has the q -> 



behavior of the exact F x . 



? EXX 



-pExx.ao i„2 -pExx.m /„ 

/q ^x /q 

pEXX,W /q pEXXSl 



(24) 



with F x °, F x °, F® 1 , and F x x , again containing a q- 
independent leading order term. 



IV. DISCUSSION AND CONCLUSIONS 

Now we analyze the physical meaning of H x and relate 
it with the GW approximation (GWA)-Bethe-Salpeter 
equation (BSE) approachJHj which represents at the mo- 
ment the most successful first-principles computational 
scheme of electronic excitations in solids. We start 
by rewriting Eq. @ as \ = (1 - Xo^xc) _1 Xo- By 
first expanding (1 — Xo-^xc) m a power series into 
l+Xo^xc+Xo^xcXo-F 1 xc + ---, next taking only the first two 
leading terms of this expansion, and finally neglecting 
correlation contributions, (1 — xoFxc)^ 1 ~ 1 + XoF x xx , 
we obtain X ~ Xo + XoF x xx xo- Thus, identifying 
XoF x xx Xo as H x [See Eq. (Q)], we can interpret H x 
as the first order correction to xo m X, 



X(q;w) « Xo(q;^) + H x (q;u). 



(25) 



Indeed, H x has been recently shown in the many-body 
diagrammatic language as the first-order self-energy .and 
vertex corrections to the irreducible polarizability x.E3 

The expression of the full H x is admittedly quite com- 
plicated. However, we point out that a simplified picture 
of the important underlying physical processes within 
TDEXX can be extracted by noting that only the first 
term of H x (H x ~ r ) and the first and third terms of H x 
(H x ~ r ) are dominant contributions at resonant w. This 
is schematically depicted in Fig. [I]. 

Note that the above situation is similar to the one 
that occurs in the solution of the BSE where the Tamm- 
Dancoff approximation is invokednJ In fact, with the 
EXX kernel, we can easily make a connection between 
the TDDFT and the GWA-BSE approach. Consider cal- 
culation, of the full response functional or excitation en- 
ergiesEllN with TDEXX and GWA-BSE. Replacing the 
bare Coulomb interaction with the screened Coulomb in- 
teraction, the resonant terms of H x effectively shift the 
EXX eigenvalue spectrum toward that of the GWA, while 
the resonant terms in H x are the counterparts of those 
occurring in the BSE in the Tamm-Dancoff approxima- 
tion. 

In the above comparison, it is interesting to observe 
that while the GWA and BSE have a clear hierarchy as 
the theory of independent quasiparticle excitationsc3 and 
electron-hole excitations, terms related to both excita- 
tion effects appear within the DFT formulation at the 
time-dependent level and the distinction between onc- 
and two-particle excitations is accordingly rather arbi- 
trary. We should also mention that the mapping be- 
tween TDEXX and GWA-BSE is not exact because the 
terms do not have counterparts in the GWA-BSE. 
These differences may indicate the inherently different 
nature of TDDFT and the GWA-BSE approach. 

In summary, we derived the expression of the EXX 
kernel for insulators and showed that it has a long- 
wavelength behavior as the exact F xc unlike the LDA 
and GGA kernels. The common conception that DFT 
is not suitable for the study of electronic excitations of 
solids was mainly derived by adopting qualitatively in- 
correct LDA and GGA potentials and kernels and the 
difficulty of going beyond them. Coupled with the al- 
ready available EXX potential, we expect the numerical 
realization of the EXX kernel will open up a new window 
of opportunity for the, first-principles study of electronic 
excitations in solids. c£l 
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FIG. 1: Schematic description of the resonant contributions to _ff x (H^ r and '). Arrow 1 and 2 represent (sk + 
q| e i(q+G) r| a j^ anc j (ak|e -! ' q+G - ) ' r |sk + q). They involve the time-sequential coupling of an electron excitation from valence 
{a, b} to conduction {s,t} bands (hole — > electron pair) and a relaxation from conduction to valence bands (electron — > hole 
pair) . 



